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Abstract 

We calculate the inelastic lifetime of an electron quasiparticle due to Coulomb interactions in an 
electron liquid at low (or zero) temperature in two and three spatial dimensions. The contribution 

of "exchange" processes is calculated analytically and is shown to be non-negligible even in the 

> 
ON ' high-density limit in two dimensions. Exchange effects must therefore be taken into account in a 
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quantitative comparison between theory and experiment. The derivation in the two-dimensional 
case is presented in detail in order to clarify the origin of the disagreements that exist among the 
results of previous calculations, even the ones that only took into account "direct" processes. 

PACS numbers: 71.10.Ay, 71.10.-w, 72.10.-d 



I. INTRODUCTION 

The calculation of the inelastic scattering lifetime of an excited quasiparticle in an electron 
liquid, due to Coulomb interactions, is a fundamental problem in quantum many-body 
theory. According to the Landau theory of Fermi liquids |lj] the inverse lifetime of an 
electron quasiparticle of energy £ p (relative to the Fermi energy Ep) at temperature T in a 
three-dimensional (3D) electron liquid should scale as 

n\{%) 2 kB T« ip «E F m w 

Te { (tIf) (" « k " T « e f • 

where ks is the Boltzmann constant. In a two-dimensional (2D) electron liquid the above 
dependencies are modified as follows 2]: 
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Besides its obvious importance for the foundations of the Landau theory of Fermi liq- 
uids pfl, the inelastic lifetime also plays a key role in our understanding of certain transport 
phenomena, such as weak localization in disordered metals. In this case, the distance an 
electron diffuses during its inelastic lifetime provides the natural upper cutoff for the scaling 
of the conductance, and thus determines the low-temperature behavior of the latter J3J, |4J, |5| . 

During the past decade some newly developed experimental techniques, combined with 
the ability to produce high-purity 2D electron liquids in semiconductor quantum wells have 
enabled experimentalists to attempt for the first time a direct determination of the intrinsic 
quasiparticle lifetime, i.e., the lifetime^that^arises purely jrom Coulomb interactions in a 

for example, the quasiparticle 
lifetime was extracted directly from the width of the electronic spectral function obtained 
from a measurement of the tunneling conductance between two quantum wells. In the case of 
large wells separation, like the ones (175A ~ 340A) studied in Ref. |8|, the couplings between 
electrons in different well are weak and can be ignored. For such weakly-coupled wells, the 
lifetime is principally due to interactions among electrons in 2D, while the contribution of 
the impurities is relatively small. 

In spite of these wonderful advances, a quantitative comparison between theory and 
experiment remains very difficult. There are several reasons for this to be so. First of 
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all, the 2D samples studied in the experiments are not yet sufficiently "ideal", namely 
disorder and finite width effects still play a non-negligible role: as a result, the measured 
lifetimes are typically found to be considerably shorter than the theoretically calculated ones. 



Secondly, the electronic density in these systems fa 
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This paper is devoted to a critical analysis of the last question, i.e., specifically, we 
calculate analytically the constants of proportionality in the relations (£Q) and (J2J in the 
weak coupling regime, and try to clear up the differences that exist among the results of 
different published calculations. One particular aspect of the confusion is the widespread 
belief that the Fermi golden rule calculation of the lifetime, based on the RPA screened 
interaction, is exact in the high density/weak coupling limit. In fact, this is only true in 3D, 
but not in 2D. To our knowledge, this fact was first recognized by Reizer and Wilkins [20], 
who introduced what they called "non golden rule processes", i.e., exchange processes in 
which the quasiparticle is replaced in the final state by one of the particles of the liquid. 
In point of truth, these processes are still described by the Fermi golden rule, provided 
one recognizes that the initial and final states are Slater determinants, rather than single 
plane wave states. In three dimensions, such exchange contributions to the lifetime were 
calculated (nomericaHy) in Ref, QQ, but the y are easily show, to beconae irrelevant in 
the high-density limit. In 2D, by contrast, the exchange contribution remains of the same 
order as the direct contribution even in the high density limit. Reizer and Wilkins found the 
exchange contribution to reduce to ~ of the direct one (with the opposite sign) in the high- 
density limit, while we find it here to be only | of the direct contribution in the same limit. 
More generally, we give an analytical evaluation of both the "direct" and the "exchange" 
contributions vs density, for boh ksT <C £ p and £ p <C ksT. 

The rest of this paper is organized as follows. In Sect. II, we provide the general formulas 
for — including exchange processes. We then devote Sect. Ill to the analytical calculation 
of — in 3D; and Sect. IV to the same calculation in 2D. The 2D calculation is presented 
in greater detail in order to explain the origin of the disagreements among the results of 
previous calculations. We explain the reason for the much stronger impact of exchange on 
the lifetime in 2D than in 3D at high density. Section V presents a comparison between 
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FIG. 1: (a) A typical scattering process between electrons of the same spin orientation near the 
Fermi surface has contributions from both a "direct" (solid line) and an "exchange" (dotted line) 
term; (b) A special class of low momentum transfer processes gives the leading-order contribution 
to the scattering amplitude in 2D at high density. 

the present theory and the experimental data of Ref. ([7() and summarizes the "state of the 
art" . 

II. GENERAL FORMULAS 



We consider an excited quasiparticle with momentum p and spin a. Its inverse inelastic 
lifetime due to the electron-electron interaction is a sum of two terms, corresponding to the 
contributions from the "direct" and "exchange" processes, respectively, 
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where £ p = 2— — \x is the free-particle energy measured from the chemical potential /i. We 
use 'D' to denote the "direct" term, and 'EX' the "exchange" term. 
Making use of the Fermi golden rule, we get |l[ , 
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and 
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where W(q) is the effective interaction between two quasiparticles, n kcr = e<3g * the Fermi- 
Dirac distribution function at temperature (3 = ^^, and we have set h— 1. The ^-functions 
ensure the conservation of the energy in the collisions. Obviously, from Eqs. (jlj and (JHJ), 
one can see that the contribution from the "exchange" process tends to cancel that from the 
"direct" process. 

As can be seen from Eq. (Jlj). there are two types of collisions contributing to the "direct" 
term, the collisions with same-spin electrons (a' = o), and those with opposite-spin electrons 
(a' = —tr). We denote the former l/r CT0 ., and the latter 1/t ct5 ., where a = —a. It can be 
easily shown that 

>-^xV (6) 
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In the paramagnetic state, one evidently has 
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The effective interaction VT(q) between quasiparticles is short-ranged compared to the 
bare Coulomb potential due to the screening effects from the remaining electrons. Such 
screening effects are normally characterized by a screening wave vector k s . Following this 
practice we approximate 
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and kp and a® are the Fermi wave vector and the Bohr radius, respectively. At very low 
density, the screening wave vector becomes much larger than the Fermi wave vector. It can 
be shown that, in this limit, 

(11) 
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or, in other words, by using Eq. (jZJ), 
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Eq. (jHj) and the low density limit result of Eq. (J12J) are exact results, which, to the best of 
our knowledge, were not explicitly established before. 

In what follows we will only consider the case of the paramagnetic electron liquid, which 
allows us to trivially dispose of the spin indices. Furthermore, by making the change of 
variable k — > k + q in the momentum summation in Eq. (jSJ), and correspondingly, k — > — k 
in Eq. (jlj), we rewrite Eqs. (JH) and © as 

— = 27T Y^ Yl W% (O^P+q^k+q 
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By using the identity, 
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where xo{<1,^) is the Lindhard function (i.e., the density- density response function of the 
non interacting electron gas), we rewrite I/t^ in Eq. (|13|) as 
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In obtaining Eq. (fTBj). we have also used the fact that 

n p+cl 5{u -£ p + £ p+q ) = 1 + e/ ^_g p / (^ - £p + ^p+q) 



Similarly, one has 
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The fact that 1/t^ and l/r^ EX ' depend only of the magnitude of p allows us to average 
over the unit vector of p = - on the right hand side of Eqs. (|16|) and (|18j) . To this end, we 
define 

n ± (q) ee ±™ - £-, (19) 
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and use the fact that 
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where 8(x) = 1 for x > 0, and 9{x) = for x < 0, and 
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Therefore l/r^ D ^ can be rewritten as 
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We note that this equation is not restricted to the regime of ksT <C Ep, but holds for 
arbitrary temperature. 

In this paper, we are only interested in the case that ksT <C Ep, and therefore the Fermi 
energy Ep is always well defined and Ep ~ ji. To perform the average over p in Eq. (|18J). 



we use the fact that, for A^T, £ p <C Ep, the contribution to l/r^ EX ^ only arises from the 
region in which £k, \u>\ <C Ep. Furthermore, the first (^-function in Eq. (fT%j) fixes the angle 
between k and q to be such as to satisfy the condition £ k — £ q +k = cu ~ 0. With this in 
mind, one obtains 

= $(p, q)6(n + (q) - uj)6{lo - Q_(q)), (23) 

where 
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A detailed derivation of this key result is presented in the Appendix. Thus finally 
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III. THE INVERSE LIFETIME IN 3D 

The theory of the electron inelastic lifetime in 3D is rather well established jjj, lHj at zero 
temperature. However, no analytical expression including the exchange has been presented 
so far, even though Kleinman ll| , and later Penn [12j , have reported numerical calculations 
of the exchange contribution. This deficiency is remedied in the present section. Our 
calculation is done at nonzero temperature, with zero temperature as a special case. 

In 3D, Eq. (g2J) becomes 
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We are interested in the case that ksT, £ p -C Ep. Therefore we only need consider the region 
of to <C Ep, in which, 
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Substituting Eq. (|27jl into 
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The integrations over q and u can be carried through, and one obtains 
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where A = ^. 

Next we move to evaluate the contribution from the exchange process. In 3D, Eq. (|25|) 
becomes 
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By using Eq. (J27|) . one has 
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After carrying out the integrations, one obtains the final result, 
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We plot the ratio of 1/t^ ex > to \/t^ d ' vs the Wigner-Seitz radius r s in Fig. 2. Notice 
that at very high density, \l/r^ EX >\ <C |l/r ( - D - ) |, and the direct-process-only theory is then 
relatively good. On the other hand, at low density, 1/t^ ex ^ = —l/2r^ D \ which agrees 
with the general conclusion of Eq. (|12|1. The contribution from exchange processes therefore 
cannot be ignored in most density range. 

In the limiting case of small excitation energy, £ p <C k B T <C E F , Eq. (|29j) reduces to 
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FIG. 2: The ratio of 1/t^ ex ^ over 1/t^ via r s in 3D. 



In the opposite of very low temperature, k B T <C £ p <C E F , one has 
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In the high density limit, A —>■ oo, Eq. (|3~4"j) becomes 
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a result obtained earlier by Quinn and Ferrell [9J 
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IV. THE INVERSE LIFETIME IN 2D 

As mentioned in the introduction, there is still some disagreement among the results of 
previous calculations of l/r e in 2D. The main purpose of this section is to exactly evaluate 
the prefactors of 1/r^ and l/r^ EX ' in 2D, and at the same time attempt to clarify the 
origin of those disagreements. We present our derivations in the two different regimes of 
k B T <C £ p <C Ep and ^ C fc^T < E F separately. For greater clarity, we also show our 
derivations for the "direct" and "exchange" contributions in separate subsections. 
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FIG. 3: The three regions of integration over q, at given uj, in Eq. 1)36(1 are labelled as I, II and III 
respectively. Only region II contributes to the leading order in 2D. 

A. ksT <C £ p : Direct process 



In 2D, for k B T < £ p < £ F , Eq 
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becomes 
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The three regions of integration, at a given value of uj are shown in Fig. 3. It can be 

shown ([22]) that the first and the third terms in the square bracket make contributions of 

the order of 0(C, 2 ), but not to the leading order of 0(C, 2 ln^ p ), which arises only from the 

second term. Hereafter we therefore focus only on the calculation of the second term. 

The small u expression for ^mxoiq^) in 2D is 

2m 2 u 
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Therefore we obtain 
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or, to the leading order, 
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Evidently, the integration in Eq. ([4*U|) has a logarithmic divergence at both the upper and 
lower limits. To the leading order, Qx{ui) can be evaluated as 

I Tfid) 
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F F 

Substituting Eq. (j41j) into (jHHjl . and performing the integration over u>, we finally arrive at 
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where we have defined the dimensionless quantity 



h,^, (42) 
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W(g) = —W(q) . (43) 

The quantity in the square brackets of Eq. (|42|) can be expressed in terms of the Wigner-Seitz 
radius r s as follows: 

W 2 (0) + \w 2 {2k F ) = 1 + 1 (y^) 2 ■ (44) 

The fact that Eq. (|40|) also has a logarithmic contribution from the upper limit of integration 
at q ~ 2^i? was missed in almost all previous analytical calculations. This is one of the main 
reasons leading to errors in the numerical prefactor of the lifetime. The second term in the 
sq ua re bracketS of Eq . m * absent in the „o rks of Ref, (fl H Q), wMle it i 8 ove, 
appreciated in the work of Jungwirth and MacDonald |l7| . where the factor | in front of 
W 2 {2k F ) is replaced by 1. All these formulas are, of course, equivalent in the high-density 
limit (r, — » 0). 

n 

Except for the work by Reizer and Wilkins [20], all the calculations cited above in 2D 
explicitly consider only the direct process, without taking account of the exchange process, 
which we deal with in the next subsection. 
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B. ksT <C £„: Exchange process 



In 2D, for k B T < £ p < E F , Eq. 



becomes 
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Again only the second term in the square bracket contributes to the leading order. Thus, 
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To the leading order, 



Substituting Eq. 
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and performing the integration over u, we arrive at 
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-W(0)[W(0) + 2W(2fc F )]ln 
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In Fig. 4, we show the ratio of l/r^ EX ^ to 1/M D ) vs r s . The remarkable fact is that, at 
variance with the 3D case, this ratio does not vanish for r s — > 0. The reason for this 
difference can be understood as follows. In 3D a typical scattering process near the Fermi 
surface, such as the one shown in Fig. 1(a), involves two particles that are well separated 
(by a wave vector of the order of 2k f) in momentum space. The direct scattering amplitude 
for such a process is maximum when the momentum transfer q is much smaller than kp, 
and is thus typically proportional to VF(0). The exchange scattering amplitude, on the 
other hand, is of the order of W(2kp) for all values of q: hence the ratio of l/r^ EX ^ to 
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FIG. 4: The ratio of 1/t^ ex ^ over 1/t( d ) via r s in 2D. 
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goes as 



W(0)W(2k F ) 

WW '■ 



which vanishes for r s — > 0. The reason why this argument fails 



in 2D is that the logarithmic contribution to the inverse lifetime in the high density limit 
does not arise from typical scattering processes, but rather, from special ones in which the 
two colliding particles are very close in momentum space (see Fig. 1(b)): hence the direct 
and the scattering amplitude are comparable, and give similar contributions to the inverse 
lifetime. A careful analysis of the integrals involved shows that in the high density limit, 
the exchange contribution cancels 4 of the direct contribution to the inverse lifetime. This 

I I 

result is at variance with that of Ref. 20j], according to which the exchange contribution 

cancels ~ of the direct one. We find that the relation 1/t^ ex ' = —\j1r^ holds only in the 

low density limit (see Eq. (|T2"|) and Fig. 4), where the weak coupling theory is not reliable. 

Combining direct and exchange contributions in a single formula we finally find that 

(50) 
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where the quantity in the square brackets is given by 




3 r s 





4 V2(r s + V2) 2 ' 
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Thus in the high density limit the total inverse lifetime differs by a factor | from the result 



of the direct-scattering-only calculation, and by a factor | from the result of Ref, 



• ([20]). 
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C. £„ <C ksT: Direct process 



For i p <tik B T <£ E F , Eq. (J22D becomes 
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where g±(k>) are the solutions of the equation fi + (g) = u>, 

g±(u;) = [p ± vp 2 — 2mw]. (53) 

Again, only the regime of to <C -Ef contributes to 1/Vp to the accuracy of the leading order. 
Thus, by using Eq. (pTfjl , one has 
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Once again only the second term in the bracket makes contribution to the leading order, 
and we find 

1 2m 3 rco 
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duj——Qi(uj) 
shpuo 



(55) 



where Qi(u) is defined in Eq. (|4())l and evaluated in Eq. (|41|1. Therefore 



m 



-(D) 



4tt 3 4 



[2iy 2 (0) + W 2 (2k F )} 



X 



which can be further evaluated leading to 

h (nk B T) 



t(°) ~ 8nE f 



uj moo 

duJ ~^~ ln T2~' 
snpuj k F 



1. 



(56) 



W 2 (0) + -W 2 (2k F ) 



In 



2E F 



(57) 
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As in the low-temperature case, the second term in the square brackets of this equa- 
tion was missed in almost all the previous theories except the one by Jungwirth and Mac- 
Donald [17], which, however, overestimates it by a factor 2. Without the second term in 
the square bracket Eq. (1571) would agree with the expression obtained by Zheng and Das 

n n 

Sarma [18|] and by Reizer and Wilkins [20(, but it would be four times smaller than the 
result of Fukuyama and Abrahams 13| , and ^- times larger than the result of Giuliani and 
Quinn 



D. £ p <C ksT: Exchange process 

In 2D, for k B T > f p , Eq. (J23J becomes 

1 m 



-(EX) 4yr 2 



rq+(u) 

duj / dq 

-oo J —q—(w) 

m+£ p /-g+M 



+ / duo I dq 

o Jq-{u) 

1 Smxo(q,oj) 



X 



sh(3u ^4p 2 q 2 - (2muj + q 



2^2 



xW{q)q[W{0) + W{\Ak 2 F -q 2 )}. (58) 



Proceeding as in the previous section we rewrite this as 

1 m 3 f°° , u 



du——Q 2 (u), (59) 



t (ex) 2 7T 3 J_ 00 sh(3u 

where C^O^) is defined in Eq. (|37)l and evaluated in Eq. (|4*%j) . Therefore 

1 m 3 



T (EX) ~ g7r 3^ 



W(0)[W(0) + 2W(2k F )} 



x w , mui 



x / du ^T ln T2" ' ( 60 ) 

/_ 00 shfJuj k F 



which can be, to the leading order, further simplified to 



I ^ kBT \w(0)[W(0) + 2W(2k F )}ln^- . (61) 



t( ex ) 32ttE f v /l v 7 v " /J 2£ F 

The ratio of 1/t^ ex > to l/r^ D ^ is therefore found to be the same as that in the case of 
k B T <C £ p , which has been plotted in Fig. 4. 

Combination of l/r (I)) in Eq. fl57|) and l/r (£X) in Eq. flM} thus yields 

- = - ( oo B P [3^ 2 (°) + 2W 2 (2k F ) - 2W(0)W(2k F )] ln ^ . (62) 

Notice the difference between the above result and the one obtained in Ref. 

16 




FIG. 5: Momentum- and energy-conserving tunneling between two identical free-electron bands 
separated by a potential difference eV is possible only if the spectral width of the single particle 
states in each band (indicated by the shaded regions) is at least as large as eV. 

V. COMPARISON WITH EXPERIMENTAL RESULTS IN 2D 



Consider two identical 2D electron liquids in closely spaced quantum wells between which 
a small potential difference V is maintained. We expect a small tunneling current between 
the layers. However, as Fig. 5 shows, no tunneling is possible in the absence of impurities 
and electron interactions. This is because under those unrealistic assumptions both the 
energy and the momentum of the electron must be conserved during tunneling, and there 
are simply no states satisfying these conditions. 

The situation changes profoundly if electron-electron interactions are allowed. Now mo- 
mentum is still conserved (if impurity scattering and surface roughness are negligible) but 
the energy of the electron quasiparticle is no longer a well defined quantity, due to the pos- 
sibility of inelastic scattering processes involving other electrons in each quantum well. As a 
result, tunneling becomes possible in a region of voltages — T < V < T where T is the width 
at half- maximum of the plane- wave spectral function in a well pj, i.e., 

i r(£ p ,T) 



A(E,Q 



(63) 



2n(E-Q2 + ir(Z p ,T)/2¥ ■ 
From the well-known relation A(E, £ p ) = — -QmG ret (E, £ p ), where G ret (E, £ p ) is the retarded 
Green's function, one can show that the spectral width T is just the inverse of the lifetime 
of a plane wave state, which is the sum of the lifetimes of electron and hole quasiparticles 
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experiment 
claculated 



r/E F 



T/T F 



Q. 



FIG. 6: Electron relaxation rate T in 2D. Experimental data are from Ref. pj, calculated ones are 
from Eq. @2J. Here T F = E F /k B . 



in the following manner:! 



m P , T) = -—^ + -j^fry (64) 



The principle of detailed balance demands 

n(€ P ,T) __ 1 -n(£ p ,T) 



(65) 



T e T h 

where n(£ p , T) is the thermal occupation number at temperature T. If we assume £ p <C k B T 
and approximate T(£ p ,T) by r(£ p = 0, T), we see from the above equations that the half- 
width at half-maximum of the tunneling conductance peak is expressed in terms of the 
electron quasiparticle lifetime as follows: 

r = 2 ^ ■ (66) 

We can now attempt a comparison between the experimental values of T from Ref. (J7() 
and the theoretical values of J-^r . This is shown in Fig. 6. It must be kept in mind that, 
in order to perform a meaningful comparison, one must first subtract from the experimental 
data a (presumably) temperature-independent constant due to residual disorder. The value 
of this constant is determined by the condition that T tend to zero for T — > 0. Even 
after this subtraction we see that the theoretical curve lies well below the experimental 
data. Furthermore, the shortcomings in Refs. [17|, [l8| as revealed in this paper imply that 
the "excellent agreement" with experiment claimed in those papers is overly optimistic, as 
pointed out earlier by Reizer and Wilkins [20]. 

We note that all derivations presented in this paper are to the accuracy of the leading 
logarithmic term. Calculations including higher order terms might bring in a better agree- 
ment with the experimental data. However, it seems too optimistic to believe that the huge 
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difference (roughly a factor 4) between theory and experiment is totally due to such higher 
order contributions. The size of the discrepancy suggests that there might be other factors 
playing a role, like the finite width of the quasi-two-dimensional system, electron-impurity 
scattering, electron-phonon scattering, and surface roughness. While the inclusion of these 
effects may help to produce better agreement with experiments, it remains a great challenge 
for experimentalists to device the conditions that will eventually allow them to probe the 
truly intrinsic behavior of the electron liquid. 
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VII. APPENDIX 

In this appendix, we give the details of the derivation of Eq. (|23j) . To this end, we denote 
the left hand side of Eq. (}2"3*j) as A s and A 2 for 3D and 2D cases, respectively. In 3D, A 3 
can be rewritten as 

A 3 = e 2 f dff sin ff f" W L + ^^ + £) 

Jo Jo V m 2m J 

1 



p 2 + k 2 + k 2 — 2pk[cos 9 cos 9' + sin 6 sin 9' cos(0 — <p')} 

(67) 

where (9, 0) and (9 f , <fi') are the spherical angles of k and p, respectively, relative to q. 
Carrying through the <$' integration, one obtains 

A, = 2ne 2 f 1 d x s( u + ?^- + ^- 
J_ 1 \ m 2m 

1 



\J(p 2 + k 2 + k 2 + 2pkcos9x) 2 — 4(p£;sin#) 2 (l — x 2 ) 

(68) 

The integral in the above equation is trivial due to the 5 function, and it leads to 

2irme 2 9(n + (q) - u)9(u - fi_(g)) 
pq\/[p 2 + k 2 + k 2 — k ■ q] 2 — [k 2 — (k • q) 2 ] [Ap 2 — q 2 ] 

(69) 



A, 
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where we have used the fact that 2muo <C k 2 Putting in this expression the approximate 

2 

equalities k ~ p ~ kp and k • q ~ — y (which follows from the condition £ic — £ q +k — w — 
due to the first 5-function in Eq. (|T%|) ) one easily arrives at Eq. ()23j) in the 3D case. 
In 2D, A 2 can be explicitly written as 



/■27T 

A 2 = e 2 / d(f)'5(uj + pq cos </>'/m + q 2 /2m) 
Jo 

1 



p 2 + k 2 — 2pk cos((f) — (p 1 ) + k s 
or, 

/•7T 

A 2 = e 2 / d(j)'5{uo + pg/m cos 0' + q 2 /2m) 
Jo 

1 

a/p 2 + /c 2 — 2p/c cos(0 — 0') + A; s 

1 

a/p 2 + A; 2 — 2pk cos(0 + 0') + /c s _ 

Carrying out the integration over <fr' yields, 

me 2 6(n + (q) - u)6{u - fi_(g)) 



(70) 



(71) 



.4, 



\/p 2 q 2 — {muj + q 2 /2) 2 
1 



p 2 + fc 2 + k ■ q - v/[A; 2 -(k-q) 2 ][4j9 2 -g 2 ] + fc, 

1 



p 2 + P + k • q + y/[k 2 - (k- q) 2 ][4p 2 - q 2 } + k s 



(72) 



Substituting, as in the 3D case, the approximate equalities, k ~ p ~ kp and k- q ~ — y one 
finally arrives at Eq. (|2*3*j) in 2D. 
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